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Abstract 

We present a NLO calculation of prompt photon production in DIS. The 
calculation involves direct, fragmentation and resolved contributions. It is 
performed in the virtual-photon proton center-of-mass system. A comparison 
of the theoretical results with HERA data is carried out. 



1 Introduction 



The prompt photon production in DIS eP collisions is an interesting reaction 
involving, in the theoretical description provided by QCD, several perturbative and 
non perturbative quantities related to virtual and real photons. At first sight the 
subprocess associated with the prompt photon production appears particularly sim- 
ple. It is the Compton effect of a virtual photon on a quark with a real photon 
and a quark in the final state: 7* + q — > 7 + q. The theoretical description of this 
process only requires the knowledge of the quark distributions in the proton and the 
calculation of Higher Order (HO) QCD corrections, the latter opening the way to a 
quantitative comparison with experimental data. 

Actually, as always with photons, the situation is quite complex. We also have to 
consider the reaction in which the initial virtual photon fluctuates into a state made 
of collinear quarks and gluons described by the virtual photon structure function. 
In the final state also, a large-p^ quark can radiate a real collinear photon, a process 
which involves a perturbative part -the bremsstrahlung of the photon- and a non 
perturbative part, and which is described by the photon fragmentation function. 
Therefore the complete description of the prompt photon production requires the 
calculation of four classes of processes. Two classes involving the virtual photon 
structure function with the final photon either directly coupled to a quark of the hard 
subprocess (the resolved-direct process) or produced through the final fragmentation 
of a parton (the resolved- fragmented process). The two other classes involve the 
direct coupling of the virtual photon to a quark of the hard subprocess with a direct 
final photon (the direct-direct process) or with a photon fragment of a parton (the 
direct-fragmented process). All four processes corresponding to these four topologies 
having, as we shall see below, the same order of magnitude. Photons can also be 
emitted by the lepton line. We do not consider this contribution in this paper. It 
can be obtained, from ref . [TJ [2] . 
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On the experimental side, HI and ZEUS [31 H] have measured the prompt photon 
inclusive cross sections for which no QCD HO calculations exist. HI also measured 
the production of a prompt photon associated with a jet. For this latter reaction a 
HO calculation does exist [2] which concerns the direct production. The fragmen- 
tation process is also taken into account and many observables are discussed in ref. 
[2]. It is the aim of this paper to present QCD HO calculations for the four processes 
of inclusive production of prompt photons in DIS eP collisions, and to compare the 
theoretical predictions with HI and ZEUS data. We do not consider the production 
of a photon and a jet. But the approach of this paper is directly applicable to this 
latter case. 

An important point in the definition of the prompt photon cross section is that 
of the reference frame in which the large-p_i_ photon is observed. There are two 
standards at HERA, the photon-proton (center of mass) CM system (hadronic sys- 
tem) and the laboratory frame. In the first one the observed large-p_i_ particle is 
produced, at lowest order, via 2 — 2 subprocesses. This frame has been used in 
almost all large-p_i_ reactions at HERA. We shall work in such a frame in the present 
paper. It is natural in photoproduction in which the almost real initial photon is 
collinear with the initial electron. For instance the photoproduction of prompt pho- 
tons, which is the Q 2 — > limit of the DIS reaction, has been studied in this frame 
[5], El [7] . A lower cut-off in p± is necessary for perturbative QCD to be valid. The 
forward production of a tt°, or a jet, in the DIS reaction has also been studied with 
this low-pt cut-off which picks 2 — > 2 subprocesses (and their HO corrections) to 
produce large-p_i_ forward partons with the aim to test the importance of the BFKL 
dynamics [HJ [91 HOI [HI [El [131 El E5]- A detailed comparison with experimental 
results was performed and a good agreement of theoretical results with the data was 
found [ni[T2l [15]. 

The situation is different in the laboratory in which a 2 — > 2 subprocess is no more 
necessary to produce a large-p_i_ particle. The transverse momentum of the observed 
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photon may come from that of the virtual photon q± = ^Q 2 {1 — y) (y being the 
inelasticity) through the basic DIS subprocess 7*+ quark — > quark, followed by the 
emission of the photon from the final quark. The authors of references [16] have 
proposed to use this frame as a way to attain the quark into photon fragmentation 
function, and have done a detailed study of this reaction at order 0(a 3 ). 

Therefore the choice of a frame amounts to emphasize a given subprocess and the 
related non perturbative quantities. To be short one could say that the laboratory 
frame emphasizes the fragmentation part, whereas the 7*-P frame emphasizes the 
virtual photon structure function and the fragmentation part. Unfortunately only 
isolated photons have been observed at HERA, corresponding to a strong suppression 
of contributions involving fragmentation functions. 

Of course it is not necessary to fix the transverse momentum and the rapidity 
of the observed photon in the 7*-P frame. This can be done in the laboratory. We 
only have to check that the subprocess in the 7*-P frame involves a large p± scale. 
For instance HI, when observing forward (in the laboratory) and large-p_i_ (i n the 
laboratory) ir°, requires a cut on the transverse momentum of the 7i~° in the hadronic 
frame [8]. We shall follow this procedure in this paper. 

In the next section we give details on the HO calculations in the hadronic frame. 
In section 3 we present results for the inclusive cross section as a function of E]_, y 1 
and Q 2 . We discuss the effect of isolation (isolated prompt photons are measured 
by HI and ZEUS). Then we examine, in section 4, the possibility to relate our 7*-P 
results with the HI and ZEUS results obtained with no p± cut in the hadronic frame. 
We will identify laboratory phase space regions in which this cut is not necessary to 
insure a large-p_i_ i n the 7*-P frame. Section 5 is a conclusion. 
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2 Technical details 



The programs containing a fragmentation function, which describe the resolved- 
fragmented and the direct-fragmented processes, can be immediately obtained from 
the programs of tt° production in DIS eP reactions [131 [15]. The only change is that 
of the fragmentation functions; we now use the Bourhis-Fontannaz-Guillet (set II) 
[IT] fragmentation functions. 

The programs describing the direct-direct and resolved-direct processes can be 
obtained from the preceeding ones by selecting the subprocesses with a final gluon 
and changing the charges and colour factors. Thus we obtain the cross sections 
corresponding to the emission of a photon. 

In all these programs the HO corrections associated with a 2 — > 3 subprocesses 
are calculated in dividing the phase space of the third final parton, which can be 
soft, in three parts: a small cylinder around the initial momenta (in the hadronic 
frame in which the virtual photon and the proton one collinear), a small cone around 
the two final hard partons, and the rest of the phase space. The cylinder allows us 
to treat the collinear and soft singularities associated with the initial partons. The 
cones allow us to calculate the collinear and soft singularities associated with the 
final partons. In the rest of the phase space the 2—7-3 cross sections have no 
singularities and the integration is performed by a Monte Carlo method [18] . 




This approach is described in details in references [T5] . It is at the root of all 
the programs of the PHOX- FAMILY [20] . What is peculiar to DIS eP reactions is 
the presence of a virtual photon structure functions. Several problems are raised by 
these functions ; for instance the implementation of the MS factorization scheme, 
their NLO evolution, their parametrizations. A detailed discussion of all these points 
has been given in reference [TJ]. Here one can only keep in mind that the available 
NLO parametrization are given in the MS factorization scheme. 

As only isolated photons are observed in the HERA experiments, we have im- 
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plemented an isolation criteria. We will use a cone criterion in the hadronic frame 
requiring no more than eEj_ hadronic transverse energy in a cone of radius R 1 sur- 
rounding the photon, with e = .1 and W = 1 (a value e = .111 will be used when 
comparing theoretical results with HI and ZEUS dateo). 

The photon transverse energy can be fixed in the hadronic frame, in which we 
require a minimum value of the latter for a perturbative approach to be valid. This 
cut also eliminates contributions from the collinear subprocess q + 7* — > q followed 
by the bremsstrahlung of a final photon. The photon kinematics can also be fixed 
in the laboratory, as it is done by HI and ZEUS. In this case, after a boost to 
the hadronic frame, the photon transverse energy may be very small. We suppress 
the contributions corresponding to this configuration by requiring again a minimum 
transverse energy in the hadronic frame. In the next section we present results 
obtained "in the laboratory frame" in order to be close to the kinematics used by 
the HERA experiments, but with a cut in the hadronic frame. 

The calculations of the photon-jet cross sections could proceed in a similar way. 
We just have to introduce a jet algorithm in the routines calculating the 2 — > 3 
subprocesses, specifying how two final partons are combined to form a jet. The 
interest of the photon- jet cross sections is due to the simultaneous measurement 
of the photon p 7 and jet pj four momenta, which allows us to define the invariant 
mass m^j = 2p 7 ■ pj. Whatever the frame in which these momenta are measured, 
perturbative QCD calculations are valid if m^j is large enough and we do not need 
any more a p_i_ 7 cut in the hadronic frame [2] . 

3 Cross sections in the hadronic frame 

In this section we present the cross sections corresponding to the four topologies 

described in the introduction, and we consider non isolated and isolated photons. 

1 The experimental cuts require that the photon carries at least 90% of the total energy of the 
jet of which it forms a part corresponding to e = .111 with our conventions. 
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We adopt kinematical parameters close to those of the HERA experiments to 
make later comparisons with HI and ZEUS easier. The beam energies of the proton 
and lepton are respectively 920 GeV and 27.6 Gev leading to \J S ep = 318.7 GeV. 
The inelasticity y = is taken in the range .1 < y < .7 and the photon virtuality 
Q 2 in the range 4 GeV 2 < Q 2 < 100 GeV 2 . (P, I and q are the four momenta of the 
proton, the initial lepton and the virtual photon). The cross section involving the 
exchange of the Z-boson are neglected. 

The rapidity of the photon in the laboratory frame is 1.8 > y 7 > —1.2 and 
the transverse energy 10 GeV > E± y > 3 GeV. In this section we do not consider 
the other cuts put by HI and ZEUS on the scattering angle of the lepton, on the 
momentum of the outgoing lepton and on the invariant mass (P + q — p 7 ) 2 = W 2 
that we shall introduce in section 4 when comparing with data. Finally, in order to 
stay in a perturbative regime, we require the photon to have, in the hadronic frame, 
a minimum transverse momentum p* ±1 > 2.5 GeV. 

We use the CTEQ6M distribution functions [23] , the parton distributions in the 
virtual photon of ref. [H] and the BFG photon fragmentation functions (set II) [1 
We work with Nj = 4 flavors. 
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Figure 1: The direct-direct cross section da/dy^ 
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Fig. 1 to Fig. 4 present the cross sections da/dy 1 corresponding to the four 
topologies discussed in the introduction. We must keep in mind that these results are 
factorization scheme dependent. They depend on the initial factorization, final fac- 
torization and renormalization scales. Only the sum of the four partial cross sections 



has a physical meaning. The scales used in this section are M = CJQ 2 + (j>j_ 7 ) 2 for 



the proton factorization scale, M 7 = JQ 2 + (C 7 j?j_ 7 ) 2 for the virtual photon factor- 



ization scale, Mp = CfJQ 2 + (pl 7 ) 2 for the fragmentation factorization scale and 



/i = C^Q 2 + (pl 7 ) 2 for the renormalization scale, with C = C 7 = Cp = 0^ = 1. 
All our calculations are performed in the MS factorization scheme. In particular the 
HO terms with a direct final photon are obtained from one-loop calculations from 
which final collinear singularities are subtracted with the MS convention. 



direct-fragmented 



Born 
Total 

Total isolated 



Figure 2: The direct-fragmented cross section da/dy^ 

Let us first consider the production of inclusive non isolated photons and let us 
start with the Born terms (Figs. 1-4). As expected the largest contribution comes 
from the direct-direct reaction. However the other reactions are not negligible. 
The direct-fragmented contribution is a factor 2 below the direct-direct one, and 
each of the resolved contributions another factor two below so that the direct-direct 
contribution represents only half of the full cross section. 
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Figure 3: The resolved-direct cross section da/dy 1 

The HO corrections are important. We display in the figures the Total = Born + 
HO cross section for each class. The HO corrections almost double the Born terms 
for the components involving fragmentation. This is due to the large number of 
subprocesses participating in the cross sections. For the contributions with a final 
direct photon, the pattern is different. For instance in the direct-direct case the HO 
corrections are very negative for y 7 smaller than .5. This produces an important 
correction of the cross section especially in the negative y 7 range. The peculiarities 
of the direct-direct reaction is the presence of only one subprocess at the Born level, 
namely 7* + q — > 7 + q, and the presence of strong kinematical constraints (parton 
distributions in the photon and fragmentation functions replaced by delta-functions). 

When the isolation is switched on, we obtain an interesting pattern of Born 
terms and HO corrections. As expected we have a strong decrease of the Born 
sections involving fragmentation functions of about a factor six, or even larger in 
the forward direction. Therefore these Born cross sections are almost negligible 
compared to the direct Born cross section on which the isolation has no effect. For 
the HO order corrections the decreases are less pronounced compared to the Born 
terms and they even increase in the cases of final direct photons; for instance by a 
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Figure 4: The resolved- fragmented cross section da/dy^ 

factor 2 for the resolved-direct case. This increase of the isolated total cross section, 
for the case where the photon is directly produced in the final state (Figs. 1 and 
3), has already been observed and discussed in ref. [24J which studies the prompt 
photon production in hadron-hadron collisions. It is due to the fact that some 
HO final collinear contributions are subtracted from the HO direct cross section 
to build the fragmentation functions. The remaining collinear terms are negative 
and the subtracted HO direct contribution increases when these negative terms are 
cut by the isolation. Of course this effect is factorization scheme dependent and 
demonstrates once more that only the sum of the cross sections has a physical 
meaning (the total cross section must decrease when the isolation is switched on). 

In conclusion, we observe that the resolved-direct isolated cross section is larger 
than the non isolated one (Fig. 3). This is also true for the direct-direct cross section 
with a very different behavior in rapidity (Fig. 1). The isolated resolved-fragmented 
and the direct-fragmented cross sections are small compared to the isolated resolved- 
direct and direct-direct contributions, but not totally negligible. 

In Figs. 5 and 6 we finally show the contributions of the resolved and the direct 
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Figure 5: Resolved and direct isolated cross sections da/dy 1 . 

processes for the cross section da/dy 1 and da/dE^ (isolated case). These figures 
illustrate the non negligible contributions of the resolved components specially for 
negative values of the photon rapidity. 

4 Comparison with experimental results 



Unlike what was done until now for the production of large p± hadrons or jets, the 
HI and ZEUS experiments do not put, in the case of prompt photons, kinematical 
cuts in the hadronic system and it is easy to check that their laboratory cuts do not 
protect the hadronic system p* ±1 against small numerical values. Thus the possibility 
of performing a perturbative calculation is not ensured. The values of Q 2 and y being 
fixed by the observation of the outgoing lepton, we have to calculate the expression 

f f dx l^t 

J dQ 2 dy J —G{x) g 4 ^ ^ (4) (p + q-P*- Ps) dp± A dy 4 dp ±3 dy 3 , (4.1) 

the definition of the momenta being given in Fig. 7a corresponding to the laboratory 
frame. 
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E Ty (GeV) 

Figure 6: Resolved and direct isolated cross sections da/dE±^. 



The delta-function can be written in terms of the light cone variables, transverse 
momenta and rapidities. 



25 (2) (q ± - p i± - p 31 _) 5 (xP (+) -^~ P±i - p ±3 e y * 



S {2) (q± - P4± - P3±) S\x- 



Q 2 y p^ + P ±3 e 



2/3 



2P-q 



P 



(+) 



5 y- 



p± 4 e m + p± 3 e 

f.(-) 



"2/3 



(4.2) 



with 



P( + ) =P n + P- 2Pr 



-) 



2E 



and =xP (+) . 

From ( 14. 2 j) we obtain the following constraints 

P_L4 



,</4 



Rj_4 = ?1 - RL3 • 



and 



(4.3) 
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We see that we can ensure a minimum value of s — 0(p±sp±4), the subprocess 



center-of-mass energy squared, with the requirement y Q 2 (l — y) = \q±\ < \p±_a\- 

To discuss the constraints on y 4 let us consider the extreme case where p 3 is 
parallel to p 4 (s — 0). With the definition q± = p^ 3 + p ±i and with y 3 = y 4 we get 
from the constraints (14.21) 

Q 2 



X = X B . 



yS, 



ep 



e y* = v _ v zL u a] 

The first result is the standard constraint associated with the ^ + 7* — > q subprocess; 
the second result shows that 2/4 cannot be large if Q 2 is small and y bounded from 
below. This can be rewritten as q± = e m 2£°y. With the HI cut y m %n — -05 we 
obtain p± m i n = 2.76 e m , which shows that the collinear configuration (s = 0) does 
not contribute to the cross section in the small p±_ and large y 4 domain. 




Figure 7: a) direct-direct kinematics in the LAB frame, b) idem for the resolved- 
direct subprocess. 



4.1 Comparison with HI data 

To determine the exact phase space domain in which our calculation is valid, 
we explore, for large value of ?/ 4 or for small values of Q 2 , the sensitivity of the cross 
section to p* ±1 - cuts. In this study we switch to the exact HI cuts which include the 
conditions E' e > 10 GeV, 153° < 6 e < 177° and W x > 50 GeV. Now the values of the 
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inelasticity y are bounded below by .05 < y, and we have 4 < Q 2 < 150 GeV 2 . We 
also implement the isolation criterion in the laboratory frame by requiring a limited 
hadronic energy E\ in a cone of radius W around the photon. This isolation is 
different from that defined in the hadronic frame since parton 3 (fig. 7) can now 
be inside the isolation cone because of the transverse boost provided by the virtual 
photon. We use the isolation parameters e = .111 and FC — 1. We perform our 
exploration of the "safe" domains by studying the direct-direct cross section. A 
typical result that we obtain for the small Q 2 domain 4 < Q 2 < 10 Gev 2 is given in 
table 1. All other HI parameters and cuts being as given above. 



Pl 7 -cut 


Born 


NLO 


2.5 GeV 


6.60 


7.14 ± .05 


1.5 GeV 


7.50 


7.64 ± .05 


.5 GeV 


7.60 


7.62 ± .05 



Table 1: Integrated direct-direct cross section in the range 4 < Q 2 < 10 GeV 2 in 
picobarns. 

We notice, as expected, a very small dependence on p*^- cut. Exploration of the 
range 10 < Q 2 < 20 GeV 2 leads to a much pronounced dependence on p* ±ry -cut, the 
cross section (Born term) varying from 4.82 to 7.2 pb. Interesting results are also 
given by HI for the range 4 < Q 2 < 40 GeV 2 , the laboratory transverse momentum 
q± = ^Q 2 {1 — y) < being not too large compared to E±^. We explore for this Q 2 
range the large y 1 range 1.4 < y 7 < 1.8 and find for the cross section a smaller 
variation then the one shown in table 1. 

To conclude we find small variations of the cross section with p^_ 7 -cut in the 
small Q 2 -range or in the large y 7 range. This incites us to compare our predictions 
with HI data in the above domains. The HI collaboration provides us with data in 
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the domain A: 4 < Q 2 < 10 Gev 2 , 3 < E ±1 < 10 GeV, -1.2 < y 7 < 1.8, in the 
domain B : 4 < Q 2 < 40 GeV 2 , 6 < £ ±7 < 10 GeV, -1.2 < y 7 < 1.8, and in the 
domain 4 < Q 2 < 40 GeV 2 , 3 < E ±1 < 10 GeV, 1.4 < y 7 < 1.8. In these three 
domains the variations with p* ±1 -cut of the cross section is small, especially domains 
B and C. In the first domain the stability of the resolved-direct (rd) cross section is 
not as good as than that of the direct-direct (dd) contribution of table I. It varies 
for 4.45 to 5.91 for p* ±1 -cut ranging from 2.5 GeV to .5 GeV. 

Our predictions are compared with data in table 3 and details of the contributions 
are given in table 2. In these three domains we set p^_ 7 -cut = 1.0 GeV. 



Domain 


dd contribution 
Born NLO 


rd contribution 
Born NLO 


df contribution 
Born NLO 


A 


1.26 


1.27 


.59 


1.18 


.0785 


.192 pb/GeV 2 


B 


.74 


.67 


.32 


.56 


.0475 


.098 pb/GeV 


C 


3.72 


4.66 


1.02 


1.88 


.153 


.370 pb 



Table 2: Details of the dd, rd and df contributions. 



Domain 


Data 


fhad 


dd + rd + df NLO 


A 


2.48 ± .21 


+.34 
-.41 


.87 


2.30 pb/GeV 2 


B 


1.78 ± .25 


+.46 
-.60 


.97 


1.29 pb/GeV 


C 


4.38 ± 1.26 


+ 1.04 
-1.75 


.77 


5.32 pb 



Table 3: Comparison between data and theory. 
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In the three domainal we find a good agreement between theory and experiment. 
Note that the last column of table 3 has been obtained by multiplying the predic- 
tions of table 2 by the hadronic correction factors fhad given by HI [3|. The small 
resolved-fragmented contribution has not been taken into account; we estimate this 
contribution to be approximatively equal to 6 % of the total cross sections. Note 
also that the contributions due to the emission of photons by the lepton line is not 
taken into account. The HI collaboration has estimated this contribution [3]. It is 
almost negligible (~ 5 %) in domain A, it represents about 13 % of the theoretical 
cross section (table 3) in domain B, and is totally negligible in domain C. 

4.2 Comparison with ZEUS data 

The ZEUS collaboration measured the isolated photon cross section at the en- 
ergy ^fs7 p = 318 GeV using the following cuts: E' e > 10 GeV, 139.8° < 9 e < 171.9° 
and W x > 5 GeV. The kinematical domain covered by the experiment is: 10 < 
Q 2 < 350 GeV 2 , 4 < E ±J < 15 GeV and -.7 < y 7 < .9. The data are quoted 
in various bins for each of these three variables with the other two variables inte- 
grated over the whole indicated range. From our previous dicussion on HI data it 
appeared that stability of our predictions under the p* ±1 cut-off could be achieved 
for low Q 2 < 10 Gev 2 , or large photon rapidity y 1 > 1.4 with moderate Q 2 values, 
however none of these domains can be extracted from the published ZEUS data. 
Restricting to large E± y values, e.g. > 10 GeV, we have tested that it is im- 
posssible to obtain stability in the full Q 2 range of ZEUS. A relative stability of 
our predictions is achieved, however, when one restricts the Q 2 range between 10 
and 20 GeV 2 integrating over the whole range for the other variables: we find a 
cross section (with only the d-d and d-f contributions) decreasing from .164 down 

2 The data in domains B and C have been obtained by subtraction between the 4 < Q 2 < 150 
and 40 < Q 2 < 150. Statistical errors are added in quadrature and we keep the largest systematic 
errors. 
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to .152 pb/GeV 2 , when decreasing the cut-off from 2.5 to 1. GeV (the d-d HO are 
not stable). Multiplying the result by a factor 2 to roughly take into account the 
resolved contribution, we obtain a cross section of .37 pb/GeV 2 (including a contri- 
bution of .045 pb/GeV 2 from the lepton line) compatible with the value of .414 ± 
.035 (stat.) pb/GeV 2 given by ZEUS. 

5 Conclusions 

A complete calculation of the isolated photon cross section in deep inelastic scatter- 
ing has been presented at the next to leading logarithmic order in QCD. Isolation 
can be imposed in the hadronic center of mass frame or in the laboratory frame. 
The calculation includes four classes of processes depending on whether the pho- 
ton is coupled directly to the hard process or through structure or fragmentation 
functions and it is valid if the momentum of the photon in the hadronic center of 
mass frame is large enough to prevent the 2— >1 process, q — > q + 7, to occur in 
a collinear configuration. Unfortunately the HI and ZEUS collaborations for their 
isolated prompt photon studies do not impose a transverse momentum cut-off in the 
hadronic center of masss frame but in the laboratory frame. This is unlike what was 
done for their studies on particle or jet production. The comparison between our 
calculation and the data can be performed only in a very restricted domain of the 
data: small Q 2 and/or large photon rapidity. It would be interesting if data could 
be available with a transverse momentum cut in the hadronic frame: a comparison 
with large transverse momentum n° production, will then be possible and it is an 
exercise always worth making as is done in hadronic colliders. The basic mechanisms 
of photon and 7r° production are different and so are the higher order corrections: 
a detailed comparison of theory with data in the case of photon production with 
many kinematical variables at hand (Q 2 , -Ej_ 7 , y 7 , £Bj) will be very constraining, as 
it was for 7r° production. 
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